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Abstract 

Atmospheric scattering of light emitted by an air shower not only attenuates direct flu- 
orescence light from the shower, but also contributes to the observed shower light. So far 
only direct and singly-scattered Cherenkov photons have been taken into account in routine 
analyses of the observed optical image of air showers. In this paper a Monte Carlo method of 
evaluating the contribution of multiply scattered light to the optical air shower image is pre- 
sented, as well as results of simulations and a parameterization of scattered light contribution 
to measured shower signal. 

1 Introduction 

Charged particles of an extensive air shower produce a large number of fluorescence and Cherenkov 
photons on their way through the atmosphere. Fluorescence light coming directly from the 
shower to the optical detector provides information needed for determining the number of parti- 
cles in the shower at points along its trajectory. The profile of longitudinal shower development 
(i.e. the number of particles in the shower versus depth in the atmosphere) is thus obtained. 
This is the principle of the fluorescence method of air shower detection. Photons emitted by the 
shower may get scattered on their way through the atmosphere between the shower front and 
the detector. Single scattering deflects the photons and therefore results in attenuation of the 
photon beam arriving to the detector. This effect is routinely accounted for in experimental air 
shower studies. However, photons may also undergo a series of scatterings before reaching finally 
the optical detector so that some of them may get redirected again toward the detector. Thus 
the atmospheric scattering results not only in attenuation of the direct fiuorescence signal, but 
may also contribute to the signal received by the detector. Since the path of multiply-scattered 
photons is longer than the direct shower-detector distance, these photons will be delayed with 
respect to direct fluorescence photons which arrive without scattering. The scattered light (both 
fluorescence and Cherenkov) must be regarded as a background for the direct fluorescence sig- 
nal, because its intensity relates to the history of shower development rather than to the current 
number of particles in the shower. 

In most cases, the optical image of a shower recorded by detectors consists mainly of direct 
fluorescence and singly-scattered Cherenkov photons. Cherenkov photons are emitted mostly at 
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small angles with respect to the shower axis, their paths are very close to the shower itself, and 
they accumulate as the shower develops. Thus singly scattered Cherenkov photons can make 
a large contribution to the shower image, as they arrive the detector simultaneously and from 
approximately the same direction as fluorescence photons, making discerning them difficult. 
A correction for the singly scattered Cherenkov photons is a standard part of all procedures 
of shower reconstruction based on fluorescence observations. Direct Cherenkov light makes a 
significant contribution to the signal in the fluorescence detector only when the shower lands 
close to the detector site or for very inclined showers [H [2]. Scattered light may be relatively 
strong in all geometrical configurations, especially in the late stages of shower development. 

Fluorescence photons, produced by the shower may also undergo scattering. A small part 
of the singly scattered photons get to the detector and are recorded together with the direct 
fluorescence light. Remaining photons, after traveling some distance in the atmosphere may 
scatter again, this time in direction of the detector. Though not likely, it is possible that a 
fraction of shower photons may reach the detector after a longer series of scatterings [3j- 

Only recently first attempts have been made to account for scattered fluorescence photons 
or multiple scattering of Cherenkov light in shower reconstruction procedures [1]. Few studies 
of this effect have been done so far. In [5j multiple scattering of fluorescence photons was 
investigated for vertical showers and an estimate of the contribution of the scattering effect 
to the shower image was obtained. An analytical study of Rayleigh scattering of fluorescence 
photons was made in [6]. A comprehensive study of scattering of light in the atmosphere is so 
far missing. 

A good estimate of the contribution of the scattering effect to the shower optical image is 
needed for precise shower reconstruction. The multiple scattering effect increases the amount 
of light arriving from the shower to the detector. Failure to account for this effect results in a 
systematic error in shower energy determination in the fluorescence method of detection. 

The experimental energy spectra determined by different experimental techniques (surface 
array and fluorescence detector) do not coincide The differences most probably are due 
to systematic uncertainties in shower reconstruction in both techniques. Systematic errors in 
the fluorescence technique are caused mainly by uncertainties in fluorescence yield, detector 
calibration, shower reconstruction and atmospheric effects - which include the effect of the 
multiple scattering. Reduction of these systematic errors is of primary importance to the studies 
of cosmic ray energy spectrum. 

The objective of this paper is the systematic study of the multiple scattering effect, by 
both Rayleigh and Mie scattering processes. A correction to the existing shower reconstruction 
procedures is obtained to account for the scattering effect. 

2 Method of simulation 

Multiple scattering of fluorescence and Cherenkov photons was simulated using the "Hybrid_fadc" 
program [8j, which was designed to simulate air shower development and detection. In this pro- 
gram, calculations are done in steps corresponding to a change of 0.04° in shower position on 
the sky, as seen by the detector. For chosen primary particle energy and shower geometry, in 
each step the program calculates the shower size using the Gaisser-Hillas parameterization and 
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Figure 1: Method of simulation (see text for details). 



the number of emitted fluorescence and Cherenkov photons. Based on these, the number of 
Cherenkov photons (both emitted directly towards the detector and singly-scattered towards 
it) are calculated. The number of fluorescence photons emitted towards the detector is also 
determined. The shower is assumed to have no lateral distribution, i. e. all photons are emitted 
at the shower axis. Calculations are done in 16 wavelength bins covering the range from 276 
nm to 420 nm. Initially, the default "Hybrid_fadc" settings were used for angular distribution 
of Cherenkov emission, molecular atmosphere profile, aerosol distribution and detector location. 
Next alternatives to these settings were tested in dedicated sets of simulations - their impact on 
the results is discussed later. 

In order to simulate multiple scattering of fluorescence and Cherenkov photons, modifications 
to the original program were made. Due to large number of photons involved, in calculations 
concerning photons produced in a given step it is impossible to trace them all separately, so 
that some simplifications are inevitable. The total number of photons is divided into smaller 
"packets" (typically 10 thousand packets in each step and each wavelength bin). All following 
calculations are done for each packet. 

A packet starts from the shower axis at the point corresponding to a current shower develop- 
ment step, with its direction randomly chosen, based on either an isotropic angular distribution 
(for fluorescence photons) or an exponential one (dN/dQ oc e~^^^^ /2'Ksin9, with — 4.5°) for 
Cherenkov light [1]. Assuming that all photons in a packet scatter at one point, two distances 
to scattering points, for Rayleigh and Mie scattering processes, are randomly drawn from the 
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corresponding mean free path distributions. From the two calculated points, the one closer to 
the starting point is chosen as the place where the scattering occurs. Knowing the geometry of 
the event and also the angular distribution of scattered light {da/dQ oc l + cos'^6 for Rayleigh or 
approximately da /dVt oc e~^/^^ ''° for Mie scattering^, respectively), the number of photons scat- 
tered towards the detector is found. Next, including the attenuation factor for the path toward 
the eye, the signal at the detector due to the scattered portion of the packet is calculated. With 
the information about the whole path in the atmosphere, the time of arrival of the scattered 
photons to the detector is found. In order to trace the remaining photons in the packet, it is 
assumed that they continue their flight together. For this smaller packet, a direction is randomly 
chosen and all calculations, just as for the first scattering, can be repeated several times. These 
calculations give as output the information about each packet: size of signal, arrival direction 
on the sky and time of arrival to the detector. 

It may happen that the point of scattering is placed by chance very close to the detector. 
In such a case, the large solid angle of the detector, as seen from this point causes the signal 
from this single scattering to be relatively very strong. In simulation results these cases can be 
recognized by the values of signal much larger than from typical packet fractions arriving from 
neighboring directions. These are artefacts of the simulation method and can be recognized in 
the results shown below. 

To investigate the instantaneous image of the shower, the signal from photons arriving simul- 
taneously, i. e. within one program step (corresponding to a change of 0.04° in shower position 
on the sky) is integrated. The duration of the program step (integration time) varies from about 

2 ns to about 40 ns, depending on shower geometry. These are short time intervals, compared 
to the signal integration time in a real detector. Moreover, in this paper we present ratios of the 
scattered to direct signal, so the variation of the signal integration time does not influence the 
results. For all time intervals of the program steps, the signal from direct and scattered light is 
calculated separately. 

3 Simulation results 

In the following, the term "new signal" on all plots denotes scattered (both singly and multiply) 
fluorescence plus multiply scattered Cherenkov light; the "shower" or "old signal" denotes direct 
fluorescence plus direct and singly scattered Cherenkov photons; the "total signal" - the sum of 
these signals. 

With the program described above, simulations were made for different shower configura- 
tions, namely for all combinations of: 

• energy - 10^^, 10^^ lO^^, lO^^eV; 

• core distance - 3, 7, 15, 25 km; 

• angle (see Fig. dj within the shower-detector plane (SDP) - 30, 50, 70, 90, 110, 130, 150 
degrees for vertical SDP; 

• SDP inclination - 30, 45, 60, 70 degrees with = 90°. 

This set of simulations was done first using the US Standard Atmosphere Model [9] and one 

^There is a large-angle tail of this distribution. Details of its parameterization were shown in [5J to have little 
effect on the overall scattering effect. 
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Figure 2: Example of simulation results (vertical lO^^eV shower landing 15 km from detec- 
tor). (A) Shower longitudinal profile and scattered light contribution in l°-radius circle around 
momentary shower position is shown. (B) The distribution of light on the sky (integrated in 
rings 0.1° wide) in shower maximum is shown. The vertical line marks the radius of the spot 
containing 90% of the signal. 
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Figure 3: Contribution of "new signal" to total shower signal versus the shower-detector dis- 
tance, in different altitude ranges (within the image spot containing 90% of signal). The "new 
signal" includes scattered fluorescence and multiply scattered Cherenkov photons. The points 
at low altitudes have largest contributions of the scattered light, and are grouped at distances 
corresponding to chosen core distances of simulated showers. 



value of aerosol concentration, corresponding to total horizontal attenuation length A7^=6.347 
km. It corresponds to a high aerosol concentration in the air - this allows us to investigate the 
effect of the multiple scattering in conditions, in which it is most prominent. Simulations with 
other aerosol concentrations were also done, and their results are presented below. 

An example of results from a single simulation run is shown in Fig. [2j Shown are contri- 
butions to light arriving to a ground detector: the intensity of light along the observed shower 
track versus its depth in the atmosphere (shower longitudinal profile). The "shower" curve 
includes direct fiuorescence, direct and singly scattered Cherenkov photons. In addition the 
contribution of scattered (including single scattering) fiuorescence photons is shown, as well as 
the contribution from multiple scattering of Cherenkov photons. The contribution of scattered 
fluorescence light (including single scattering) is at all stages larger than from multiple scat- 
tering of Cherenkov light. The signal from multiply scattered light is larger at later stages of 
shower development, and may finally reach a few percent of the total signal from the shower. 
To compare contributions of scattered light in various shower geometries, in Fig. [3] shown are 
contributions of scattered fiuorescence and multiply scattered Cherenkov light within the "image 
spot", i. e. a circle containing 90% of light from the air shower image [2j. 

Results from the whole set of simulations performed are presented in Figures [3] and |H Shown 
are contributions from multiply scattered light for all showers at their maxima and points 250 
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Figure 4: Contribution of singly and multiply scattered fluorescence (A) and multiply scattered 
Cherenkov (B) to total shower signal plotted versus altitude above ground (within the image 
spot containing 90% of signal). 
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g/cvc? before and 250 g/cvo? after the maximum (if above ground). In Fig. [3] the simulation 
points group around 3, 7, 15 and 25 km distance because of chosen core distances in simulations. 
The relative contribution of the scattered light appears to be strongly correlated with altitude 
above ground, rather than with distance from shower to detector. This contribution increases 
with decreasing altitude. Figure [3] suggests therefore that the altitude above ground rather than 
the shower-detector distance may be a better parameter to organize the data. The same results 
are plotted in Fig. [Has a function of altitude. The contribution of scattered fluorescence photons 
and multiply scattered Cherenkov photons are shown separately. The contributions from shower 
maxima and points earlier and later in shower development that are at the same altitude, show 
also comparable contribution of scattered light - no clear dependence on shower age is seen. 
On the plots shown are relative contributions of multiply scattered light to the shower image, 
within the 90% "image spot". When calculated in this form, the scattering contribution does 
not depend on the size of the shower image nor on its brightness. The stochastic nature of 
the scattering process, enhanced by the packet algorithm and limitation of available simulation 
time, lead to a relatively large spread of points. 

To investigate the impact of the vertical distribution of the atmosphere (the so-called molec- 
ular atmospheric profile) on the contribution of multiply scattered light, smaller sets of simu- 
lations were made. Compared were the US Standard Atmosphere Model (used in all previous 
simulations) and the models of atmosphere in January and July at the southern Pierre Auger 
Observatory in Malargiie (Argentina) [10^ [TT] . Results of the simulations show that changing the 
molecular atmospheric models has little influence on the final results, as shown in Fig. [5]A.: the 
effect of changing the atmospheric profile is smaller than the scatter of points due to variation of 
other parameters. In other words, the dependence on the atmospheric profile is a second order 
effect. 

Compared were also results of simulations using different angular distributions of Cherenkov 
photons emitted by a shower: simple, one-exponential distribution: clN/dil. oc e~^^'^'^° /lirsinO 
with a more realistic two-exponential one jT2] : 

dN/dn oc e-^^'^-^° /27rsin9,9 < 35° and dN/dQ oc e-^/'^'^-^° /2TTsine,0 > 35° 

The results are shown in Fig. [5f3. Again, no significant difference between these two sim- 
ulation sets can be seen. Both models of Cherenkov emission differ only in distribution of a 
very small fraction of photons emitted at large angles. On the other hand, the distribution of 
photons after a series of scatterings is not expected to be very sensitive to fine details of the 
original angular distribution of only a small part of the emitted photons. Thus, the effect of 
details of Cherenkov emission distributions appears to be a second order effect. 

All described above shower simulations were done for a detector placed at an altitude of 1570 
meters above sea level, which is the default altitude in the Hybrid-fade program and corresponds 
to altitude of HiRes detectors (1550 and 1593 m a.s.l.). It is also roughly the mean of the alti- 
tude range of fiuorescence detector locations at the southern Pierre Auger Observatory (located 
between 1421 and 1719 m a.s.l.). It is necessary to check if the scattering contribution changes 
with altitude of the detector, since the air density, important for calculations of scattering, 
changes exponentially with altitude. It should be noted that Rayleigh scattering changes with 
air density, that is with altitude above sea level, while aerosol concentration, and consequently 
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Figure 5: Comparison of simulation results for different atmospheric models (A) and different 
angular distributions of Cherenkov emission (B). 
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detectors above sea level. 



Mie scattering, changes independently according to local conditions and altitude above ground. 
Again, sets of simulations were made, with the detector placed 150 meters below and 150 meters 
above the default altitude, i.e. at 1420 and 1720 meters a.s.l. The results are shown in Fig. [6l 
Different simulation sets show no significant difference in the scattering contribution. It may 
not be surprising - the difference of 300 meters in altitude translates to only 3% change of air 
density (and in horizontal attenuation length for Rayleigh scattering), so no large changes in 
the scattering effect should be expected. 

The results presented above demonstrate the contribution of the multiple scattering effect 
to the optical image of the shower defined as a circle containing 90 % of the signal. In a real 
detector, the field of view is divided into pixels of a fixed size (1.5° in diameter in the Pierre 
Auger Observatory fluorescence detector [13], 1° in the HiRes detector |14j). In the shower 
reconstruction procedure, the shower signal is obtained by summing signals from pixels located 
within some angle C from a current center of the shower image (Fig. [7]). The size of the C angle 
is chosen for each shower individually, to maximize the signal-to-noise ratio [13j. On Fig. [8]A. 
shown is the dependence of the C angle on the distance from the shower to the detector, for 
showers recorded by the Auger Observatory [15j. The exponential function fit to the data is 
compared to the size of the 90% image spot in Fig. [Sf3. For the distant showers, the C angle 
is larger than the radius of the image spot. Detector properties, in particular pixellization of 
the field of view, are important here - the C angle has to be always larger than the radius of a 
pixel. On the other hand, the spot size of the optical image of a shower depends on geometry 
and decreases with distance: being comparable to C angle for nearby showers, the image size 
becomes much smaller for distant showers. 
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Figure 7: Schematic illustration of an air shower observation in a detector. The field of view of 
the detector is divided into hexagonal pixels, each observed by a photomultiplier. Light from 
air shower triggers pixels along the observed line of its propagation. The recorded signal is 
integrated within different angles (" to find the best signal-to-noise ratio. 



As demonstrated in Fig. [2j3, the angular distribution of scattered light is much wider than 
that of the direct light (i.e. wider than the "shower" curve). Therefore, the relative contribution 
of the scattered light to the recorded signal depends on the C angle. Since the goal of this work 
is to account for the scattering effect in analyses of data from real detectors, the dependence of 
the scattering contribution on ( is essential. 

For nearby showers, C is comparable to the 90% spot radius. With angle larger than 1°, for 
distant showers, the light from solid angle much larger than shower image spot is recorded in the 
detector. This means that the contribution from multiple scattering to the recorded signal in the 
detector becomes larger for distant showers than for nearby ones. The relative contribution of 
the scattered signal calculated not for the 90% spot radius (as in Fig. [3]), but for various angles 
is shown in Fig. [U] and [TUl The scattering contribution grows roughly linearly with increasing 
distance (this can be seen by comparing scales of vertical axes on Fig. [TOll . For typical C values 
of 1°-1.5° used in shower reconstruction and for distant showers, the "new signal" from scattered 
fluorescence and multiply scattered Cherenkov photons exceeds 10% of the direct shower signal 
at low altitudes. As shown on plots of Fig. [TOl the character of the dependence on distance 
and altitude is similar for different values, and the value of the scattering contribution scales 
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Figure 8: (A) Averaged values of C angle used in shower reconstruction (from Auger data |15j ) 
for different shower distances and an exponential fit to these points. (B) This C fit compared 
with the radius of image spot containing 90% of shower light. 
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Figure 9: Contribution of new signal to shower signal versus shower-detector distance within 
light collecting angle C = 1°. 



approximately linearly with On all plots of Fig. [10] presented are simulated data points 
representing maximum of shower signal, and also at 250 g/cm^ before and 250 g/cm^ after 
the shower maximum. No significant separation of data points at different stages of shower 
development (i. e. different shower age) can be observed; all points at a given altitude and 
distance have comparable values of the scattering signal, so the scattering contribution appears 
to be independent of the shower age. 

Another variable that must be taken into account is the amount of aerosols in the air. All 
previous calculations were done with one value of aerosol concentration, equivalent to horizontal 
attenuation length (at ground level) for Mie scattering of Am=9.59 km, which is the default 
value of the "Hybrid_fadc" program. For the Rayleigh scattering at the default detector altitude 
(1570 m a.s.l.) the horizontal attenuation length for the light wavelength of 361 nm is A/j=18.77 
km. Therefore, the total horizontal attenuation length (l/Ay = I/Ar + 1/Am) is A'j'= 6.347 
km. A/j and Am cannot be measured separately. Ar depends on the air density only, so it is 
well known. Am can be obtained only by measuring A^, so the total attenuation length At 
can be used as a parameter in study of aerosol concentration variations. In order to check the 
influence of amount of aerosols on the contribution of multiple scattering, a set of simulations 
was done: 

• total horizontal attenuation length At- 9.487, 11.827, 15.699 km (corresponding to 50, 30 
and 10 percent of "Hybrid_fadc" default aerosol concentration); 

• shower energy: lO^^eV; 
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different values of light collecting angle Q. The lines represent fits of Equation ([T|) for d = 4, 7, 
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• shower core distance: 3, 7, 15, 25 km; 

• angle: 30, 50, 70, 90, 110, 130, 150 degrees for vertical SDP; 

• SDP inclination: 30, 45, 60, 70 degrees with ^ = 90°. 

Examples of results from these simulations are shown on Fig. [TTl It can be observed that 
the contribution of multiple scattering varies strongly with aerosol concentration. Especially for 
low altitudes above the ground, a dependence is seen: for higher aerosol concentration (smaller 
At) the scattering contribution is higher. 

To describe fully the concentration of aerosols in the whole volume of air, not only Km must 
be known, but also the vertical distribution of aerosol particles. Generally it is assumed that the 
concentration of aerosols falls down exponentially with altitude, with a scale height of 1.2 km. 
This distribution was used in all multiple scattering simulations discussed above. In order to 
see if a change of aerosol scale height influences the final results, a set of simulations was made: 

• aerosol scale height Ha/: 0.6, 2.4 km: 

• At = 9.487 km; 

• energy lO^^eV; 

• core distance: 3, 7, 15, 25 km; 

• angle: 30, 50, 70, 90, 110, 130, 150 degrees for vertical SDP. 

Results from these simulations are shown on Fig. [T2l The data points were selected so that 
they all have comparable values of optical depth r for the line of sight from the emission point at 
the shower to the detector. The optical depth was selected because it describes both geometrical 
distance and scattering properties of the path from the air shower to the detector. On the plots 
compared are the multiple scattering contributions for different aerosol scale heights. It can be 
observed that even large differences of aerosol scale height (within a factor of 4), don't make 
a significant change - no regular dependence on the scale height is seen. The contributions of 
multiple scattering for the same optical depth and altitude show some variation (1-2%). With 
the aerosol scale height changed, the scattering contributions are still within the range of this 
variation. Therefore a conclusion can be made that variation of scale height of aerosols doesn't 
need to be separately accounted for. 



4 Parameterization of scattering contribution 

All observations made during the review of simulation results allow one to make a parameteri- 
zation of the multiple scattering contribution. This contribution is calculated as: 

new signal 
old signal 

The old signal is the sum of direct fiuorescence, direct Cherenkov and singly scattered Cherenkov 
light. The new signal is the sum of scattered (singly and multiply) fiuorescence light and 
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Figure 11: Contribution of the new signal to shower signal versus altitude above ground for 
different aerosol concentrations characterized by different values of total horizontal attenuation 
length At- The contributions are integrated within the angle = 1°. Shown are points at 
selected distances from the detector. 
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Figure 12: Comparison of the new signal for different values of aerosol scale height Hm- On the 
plots are shown points with similar values of optical depth. The contributions are integrated 
within the angle C = 1°. 
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Figure 13: Contribution of the new signal to shower signal versus altitude above ground for 
selected shower-detector distances for = 1.5°. Solid lines represent fits to function ([T]) to all 
simulation points for d = 4, 7, 15 and 25 km. Dashed lines show field of view limits of Auger 
fluorescence telescopes (2 and 30 degrees above the horizon). 



multiply scattered Cherenkov light. For one value of horizontal attenuation length A^, M may 
be parameterized as a linear function of the ( angle (in degrees) and of the shower-detector 
distance d and an exponential function of altitude above ground h (both in kilometers)|l: 

M = ACdeM-^) (1) 

A fit of this function with three independent variables {(,d,h) and two parameters {A, B) was 
made. Data from simulations for Aj'=6.347 km (at 361 nm) were used, representing shower 
image at their maxima and points 250 g/cvc? higher and lower (if above ground), in the range 
of the Q angle between 0.1° and 5° with a step of 0.1°. The resulting parameters of the fit 
are: A = 0.5830 ± 0.0011%/(iec/, B = 2.4986 ± 0.0062km, when C is in degrees, d and h in 
kilometers. 

In order to confirm that the dependence on C, and d is linear, an additional fit to a following 
function was made: 

M = AC^d^exp(-4). 

B 

The resulting values of the exponents: C = 0.986±0.007, D = 0.993±0.006 show that treating 
the scattering contribution as a linear function of C, and d is an acceptable simplification. 

^We note that two geometrical variables are needed to make the parameterization: we chose d and h 
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Fits of Eq. [T]to the simulation results are shown in Figure [TOl To better show the agreement 
with the data points, a subset of the data in smaller intervals of distance d is shown in Figure 
[T3l The fits are shown for d = 4, 7, 15 and 25 km, with data points corresponding to distances 
within ±1 km from these values. Equation [T] fairly well represents the contribution of the new 
signal to shower signal in a fluorescence detector. 

The scattering contribution strongly depends on altitude of the shower front above the 
ground, and consequently - on the elevation angle of the line of sight. If the field of view 
of a detector is limited, it may limit the range of the scattering contribution values that are 
really observed. As an example, the approximate limits of the field of view of the Pierre Auger 
Observatory fluorescence detectors (between 2° and 30° above the horizon) are marked in Fig. [13] 
by the dashed lines. This means that only points located to the right of the "elevation=2°" line 
represent contributions which are relevant for the Auger detectors. Nevertheless, the multiple 
scattering contribution to the recorded signal may exceed 10%. 

For other values of the horizontal attenuation length, similar fits can be performed. The 
results for two of the other A^^ values are shown on Fig. [TH As it can be observed. Equation 
[T] describes fairly well the scattering contribution for different aerosol concentration, especially 
for relatively large values of the contribution. However, the A and B parameters of Equation 
[T] depend on the horizontal attenuation length (see fig. [TS]) . Therefore, in order to use this 
function as parameterization of multiple scattering contribution, the A and B parameters must 
be described as functions of At- And so, as it is shown in Fig. \T5\ parameter A can be 
approximated by an exponential function, and B - by a linear function. When presenting the 
Equation [T] in this form, the following values are obtained from the fit: 

A = ai exp(-AT/a2) + «3, 
where m = 1.774 ± 0.033%, 
02 = 4.365 ± 0.064km, 

as = 0.1387 ± 0.0025%; (2) 
B = biAx + b2, 
where h = 0.1976 ± 0.0035, 
b2 = 1.402 ± 0.029km. 

With these expressions for the A and B parameters, the parameterization of the multiple 
scattering contribution by Eq. [T]and Eq. [2] is complete. This contribution is given as a function 
of the signal integration angle C; the shower-detector distance d, the altitude of the current 
shower position h and the total horizontal attenuation length At- 

For different aerosol concentrations in the air, the scattering contributions are different, even 
for the same geometrical distance and altitude of the shower front. Therefore it can be concluded 
that a parameter independent of individual geometry could be more convenient. The optical 
depth r is a parameter that characterizes the distance between two points, with respect to the 
atmospheric scattering on the way between these points. Indeed, an observation can be made: 
if the scattering contributions for different At are grouped not by their geometrical distance 
d but by optical depth r, the scattering contribution values are similar for the same altitudes. 
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Figure 14: Contribution of the new signal for different values of the total horizontal attenuation 
length. The contributions are integrated within the angle C = 1°. The lines represent fits of Eq. 
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Figure 15: Dependence of the A and B parameters (Eq. [T]) on the total horizontal attenuation 
length. 



Therefore, the function to which the data is fitted, may be the following: 

M = FCrexp(-^) (3) 

r is calculated for 361 nm wavelength. A fit to the function [3] gives the following values of the 
parameters: 

F = 3.3179 ± 0.0078%/dec/, G = 5.426 ± 0.027km. (4) 

As it is shown in Fig. [TBI a function of this form also gives a reasonably good fit to the results 
from simulations. The fitted function describes the results from simulations with accuracy of 
about 1-2 %, which is sufficient for applying this correction. Statistical fluctuations of results 
from individual simulations don't allow a more accurate parameterization. In this form the 
parameterization of multiple scattering contribution is a function of the C angle, optical depth 
r along the shower-detector line and the geometrical altitude of shower h. 

The results of our simulations are compared to results of other studies [U [6], wherever 
possible. In Ref. |5j the multiple scattering of fluorescence light for vertical showers was studied, 
assuming a uniform fluorescence light intensity along the shower path. To make the comparison, 
scattered fluorescence light only from the vertical showers from our simulations is plotted in Fig. 
[T71 as a function of quantity used in [5| : t * a * y/d, where r is the optical depth between the 
shower and the detector, a is the total scattering coefficient at the point of emission [m~^], d 
is the shower-detector distance [m]. Results for ( = 1° and C = 2° are compared. Although 
the methods of simulation in both studies are different, the results appear to be in a reasonable 
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Figure 16: Contribution of the new signal for all values of horizontal attenuation length. Points 
are grouped by their value of optical depth for the shower-detector distance. The lines are fits 
of Eq. [3l drawn for the mean value of respective r ranges. 



agreement, especially for small One can observe on Fig. [T7] that the contribution from 
scattering grows as a function of r * a * ^/R, except for the final phase of shower development. 
This can be explained by proximity to the ground level - larger part of light produced by 
shower particles is absorbed at the ground before scattering in air can occur, thus decreasing 
the scattering contribution. 

Another comparison was made of our simulations with analytical calculations of Ref. [6]. 
These analytical calculations are done for Rayleigh scattered fluorescence light, assuming a 
constant-density atmosphere. To make the comparison, our simulation program was modified - 
allowed was Rayleigh scattering only. To simulate the effect of the constant-density atmosphere, 
for comparison were chosen points low above the horizon (2°). However, one must remember that 
in a 3-dimensional simulation the photons traverse layers of atmosphere with different densities; 
also proximity to the ground level (elimination of photons at altitude zero) may influence the 
results. Although both calculations were done using completely different methods, final results 
are in a satisfying agreement (Fig. [TH]) . The analytical formula of [6J for correction is also 
consistent with our fitted function - the multiple scattering contribution scales linearly with 
angle and distance. 
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Figure 17: Comparison of fluorescence scattering from our set of simulations (data points) with 
results of [5] for ^ = 1° and Q = 2° . The lines represent the parameterization given in Ref. [5]. 
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Figure 18: Comparison of results from our simulations of Rayleigh scattering only with the 
analytical function 3.1 * k * ( (where k = d/X) from [6], for different shower-detector distances 
and wavelength bins. 
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5 Conclusions 



A Monte Carlo method for simulation of the multiple scattering of light has been developed, and 
used to determine the effect of this new contribution on the observation of air showers. Simu- 
lated is the development of air showers in a realistic atmosphere. Calculated is the contribution 
of the new signal: scattered (both singly and multiply) fluorescence photons and multiply scat- 
tered Cherenkov photons to the observed shower image. The simulation results show that the 
new contribution to the signal recorded by a detector scales linearly with the signal collecting 
angle C and the shower-detector distance (geometrical distance d or optical depth r), and falls 
exponentially with the shower front altitude above the ground; no dependence on the shower 
age is seen. 

Different distributions of molecular atmosphere and aerosols, or angular distributions of light 
emission by an air shower were used in the simulations. It can be concluded that the size of the 
new signal depends on the distance between the shower and the detector, on the altitude of the 
observed point above ground, on the signal collecting angle in the detector and on the amount 
of aerosols in the air. It has been shown that one doesn't need to account separately for the 
variations of the molecular atmosphere profile, the vertical scale height of aerosol distribution, 
Cherenkov angular distribution or different altitudes (above sea level) of the detector. 

As a result of this study, two parameterizations of the multiple scattering contribution to 
shower signal have been obtained. These parameterizations are simple functions of parameters 
characterizing the location of the air shower and the atmospheric conditions. One parameteri- 
zation (Eq. [T]with parameters given in Eq. [2|) is a function of the signal collecting angle the 
shower altitude above ground h, the shower-detector distance d, and the horizontal attenuation 
length At (at 361 nm). The other parameterization (Eq. [3] with parameters given in Eq. S]) 
allows one to calculate the new signal based on (, h and the optical depth along the shower- 
detector line of sight in the atmosphere r. These parameterizations can be easily implemented 
into existing reconstruction procedures used in analyses of air shower observations. 

Scattering of fluorescence light and multiple scattering of Cherenkov photons makes a con- 
tribution to the signal received by the fluorescence detector, that should be taken into account 
in analysis of experimental data. This contribution can reach and exceed ~ 10% of the shower 
signal. Since the multiple scattering contribution changes with altitude, applying a correction 
for this effect will slightly change the shape of shower profile (and with it the depth of shower 
maximum) . 

In the form as it is presented in this work, the scattering correction should be included into 
reconstruction procedures after geometry of air shower and signal received at the detector are 
calculated. Another correction for the lateral distribution of light in shower image [16j is also 
applied at this stage. (This correction accounts for the part of shower signal that is lost outside 
the signal collecting angle (.) We may notice that, while these two corrections have opposite 
effects on reconstruction results (taking multiple scattering into account reduces the signal from 
a shower, whereas lateral distribution correction increases it), they only rarely cancel each other: 
the effect of multiple scattering is larger for the distant showers than for the nearby ones, which is 
opposite to the lateral distribution correction. The application of both corrections is important, 
as they decrease the systematical uncertainty of final results. 

It should be considered that in future not only shower reconstructions, but also simulations 
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should include the effect of atmospheric multiple scattering. Also, in order to get consistent 
results, a number of calibration procedures may have to be amended to account for multiple 
scattering. Some methods of detector calibration and measurements of atmospheric conditions 
are based on observations of a distant light source or laser beam. Results of such measurements 
can also be influenced by the multiple scattering. 
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